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ABSTRACT 

We have tested the reliability of the red giant branch (RGB) as a metallicity indicator ac¬ 
counting for observational errors as well as the complexity of star formation histories (SFHs) 
and chemical evolution histories observed in various stellar systems. We generate model color- 
magnitude diagrams (CMDs) produced with a variety of evolutionary histories and compare the 
resultant metallicity estimates from the colors and magnitudes of RGB stars to the true input 
metallicities. We include realistic models for photometric errors and completeness in our synthetic 
GMDs. As expected, for simple stellar populations dominated by old stars, the RGB provides a 
very accurate estimate of the modular metallicity value for a population. An error in the age of 
a system targeted for this type of study may produce metallicity errors of a few tenths of a dex. 
The size of this metallicity error depends linearly on the age error, and we find this dependence 
to be stronger with more precise photometry. If the population has experienced any significant 
star formation within the last ~6 Gyr, the metallicity estimates, [M/H], derived from the RGB 
may be in error by up to ~0.5 dex. Perhaps the most important consideration for this technique 
is an accurate, independent estimate of the average age for the target stellar system, especially 
if it is probable that a signihcant fraction of the population formed less than ^6 Gyr ago. 

Subject headings: stars: abundances - stars: HertzsprungRussell and CM diagrams - galaxies: dwarf - 
galaxies: stellar content - Local Group 


1. Introduction 

The red giant branch (RGB) phase in intermediate- 
to low-mass stellar evolution corresponds to the 
phase just before the onset of helium fusion. Dur¬ 
ing the RGB phase, a star contains an electron de¬ 
generate core surrounded by a hydrogen-burning 
shell fueled by the GNO cycle ( Ghiosi|[l998 1. 

The color of the RGB on a color-magnitude dia¬ 
gram (CMD) strongly depends on the metal abun¬ 


dance of the stellar population (e.g. Da Gosta & 


Armandroff (1990)). For this reason, the RGB 
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has been widely used to obtain average metallic¬ 
ity estimates for a variety of stellar systems in the 
form of straightforward empirical calibrations be¬ 


tween magnitude, color, and metallicity (Da Gosta 


fc Armandroff|1990{ Saviane et al.|2000[|Streich et 


al. 2014). However, the morphology and position 


of the RGB also depend somewhat on the age of 
the population, especially for stellar populations 
younger than ~10 Gyr. Thus, given an estimate 
for the age of a stellar population, the metallic¬ 
ity distribution function (MDF) of the system can 
be calculated by comparing the RGB with theo¬ 


retical isochrones of the appropriate age (Saraje- 


dini fc Jablonka|2005 McConnachie et al. 20061 or 
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fiducial RGB sequences for discrete ages (]Durrell, 


Harris, & Pritchet||2001 

I- 



Previously, 

Salaris & Girardi 

(2005 

) had in 


vestigated biases in determining the tip of the 
red giant branch (TRGB) based in the galaxies 
LMC, SMG, and LGS3 using synthetic GMDs. 
That study investigated some population effects 
on metallicities derived from the synthetic GMDs 
representing these systems. In this work, we pro¬ 
vide a more thorough investigation of the de¬ 
gree to which star formation history can influ¬ 
ence the results of these types MDF analyses men¬ 
tioned above. To facilitate this, we generate model 
GMDs with a variety of star formation histories 
(SFHs) and chemical evolution histories and cal¬ 
culate the MDFs by interpolating amongst a grid 
of isochrones. In Sectionj^ we outline the methods 
used, detailing how the synthetic GMDs were cre¬ 
ated, how error and completeness profiles were ap¬ 
plied, and how the interpolations were performed. 
Results are presented in Section and we con¬ 
clude our findings in Section]^ 

2. Methods 


2.1. Synthetic CMDs 


To genera te the synthetic CMDs, we utilize 
lAC-STAI^ ( [Aparicio fc Gallart|[2004 ). This web 
browser-based software package generates syn¬ 
thetic CMDs through bilogarithmic interpolation 
in age and metallicity for a given set of stellar 
evolution libraries. lAC-STAR allows the user to 
enter unique SFHs and age-metallicity relations 
(AMRs) to synthesize arbitrarily complex stellar 
populations. It does so by taking as input up 
to 20 nodes in time with corresponding star for¬ 
mation rates (SFRs) and metallicities. The code 
then performs interpolations between these nodes 
to achieve satisfactory temporal resolution. 

The stellar properties are calculated for each 
star and then converted to absolute magnitudes 
in a range of different filter sets with a variety 
of bolometric correction libraries. We utilize the 


stellar evolution library of Girardi et al. (2000) 


along with the bolometric correction library from 
Girardi et al. (2002). We chose this combina¬ 


tion after our own exploration revealed that the 
RGBs produced from lAC-STAR with the stel¬ 


lar evolution library of Girardi et al. (20001 more 


closely matched the corresponding isochrones than 
the other stellar evolution libraries made available. 
We retrieved the isochrones to construct the in¬ 
terpolation grid using the CMD v2.1 web inter- 
fac^E] since these utilize the same bolometric cor¬ 
rection library used by lAC-STAR. For simplicity, 
we use a distance modulus of zero, zero redden¬ 
ing, a Kroupa IMF ( Kroupa||2001 ), and an age of 
13.5 Gyr for the universe in all of our models. We 
also set the mass-loss parameters for the RGB and 
AGB to f] = 0.2 since we focus on old (low-mass) 
populations. As we are interested in magnitudes 
consistent with the brighter RGB (—4 <Mi < 0) 
in the present study, we set the limit for output 
from lAG-STAR to My = 2 mag. 


2.2. Photometric Errors and Complete¬ 
ness 

In order to realistically model the observational 
characteristics of GMDs, we follow the prescrip¬ 
tion of Barker, Sarajedini, & Harris (2004). In 


that work, they tested the reliability of the TRGB 
using synthetic GMDs where they modeled the 
photometric error and completeness profiles using 
simple analytic functions. Photometric errors for 
both My and Mi are modeled with an exponential 
function: 

(t(M) = (1) 

where k and r are coefficients that describe the 
shape of the error profile. The completeness at a 
given magnitude is modeled as: 


/(M) =-arctan[a(M — Mq)] (2) 

TT 

where a is a shape parameter and Mq is the mag¬ 
nitude at which completeness falls to 0%. 

In an effort to sample two extremes in error and 
completeness, we define two error and complete¬ 
ness profiles in the following way. For the first 
case, we utilize the same error and completeness 


profiles from Barker, Sarajedini, & Harris (2004), 


which will be denoted as observational profile A. 
The relevant coefficients for this error profile are 
listed in Table This results in typical errors 
on the RGB of <0.05 mag at absolute magnitudes 
brighter than M = -2 mag and a completeness of 
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^ http: //stev.oapd. inaf. it / cgi-bin / cmd_2.1 
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Fig. 1.— Both adopted error and completeness pro¬ 
files. Top: Photometric error as a function of absolute 
magnitude, both in V and I. The dashed line corre¬ 
sponds to Profile A, while the solid line corresponds 
to Profile B. Bottom: Completeness as a function of 
Mv. Symbols are as for the top panel. 


50% at My = 1 mag. A typical CMD obtained 
using this profile contains ~2200 stars. 


The second case is aimed at representing a less 
conservative error profile similar to Hubble Space 
Telescope (HST) imaging of a Local Group dwarf 
galaxy, which we will denote as observational pro¬ 
file B. In particular, we utilized the photometry 
from Hidalgo et al. (2009) using HST imaging 
of the dwarf transition-type galaxy, Phoenix, for 
which the data are deep enough to sample the 
RGB 6 magnitudes fainter than the TRGB. From 
the photometry provided by S. Hidalgo (private 
communication; 2013), we fit functions of the form 
shown in Equation 0 and ([^ to the photomet¬ 
ric errors and completeness. We note that [Barker, | 
Sarajedini, fc Harrisj (|2004|) utilized the same error 


profile for both V and I, and we do the same here, 
as we found there to be little difference between 
them. We therefore use the fit from the I-band 
photometric errors. The coefficients for profile B 
are also listed in Table for comparison. Typical 
errors on the RGB for Profile B are ~0.01 mag. 
Both error and completeness profiles are plotted 
for comparison in Figure Finally, in order to 
realistically simulate the density of stars on the 
CMD for a system like a dwarf galaxy, we scale 
each synthetic CMD such that the number of stars 
at Mi = —3.0 ± 0.1 mag is ~90. A typical CMD 
obtained using this profile contains ^3500 stars 
brighter than My = 2 mag. 


2.3. Metallicity Calculation 


In an effort to derive metallicities from these 
synthetic CMDs in a manner similar to what 
has been done in past work, we chose to calcu¬ 
late the metallicity of each individual star on the 
RGB by interpolating among a grid of theoreti¬ 
cal isochrones. Our procedure is modeled after 


that of Sarajedini & Jablonka (2005), and we uti¬ 


lize the same Interactive Data Language (IDL) 
code to perform the interpolations. To summarize 
the procedure, the isochrones represent a three- 
dimensional grid of Mi, (V-I)o, and metal abun¬ 
dance for stars of a given age. The IDL routine tri¬ 
grid then performs the interpolation for each star 
on the grid within a specified magnitude range. To 
estimate the errors in metallicity arising from the 
presence of photometric errors, we repeated the 
interpolation to account for photometric error in 
magnitude and color for each star, and the sum of 
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Table 1 

Coefficients for different observational profiles. 


Profile 

n 

T 

a 

Mo 

A 

0.11 

0.39 

2.0 

1.0 

B 

0.02 

0.67 

1.18 

4.57 


these deviations in quadrature is assigned as the 
error in metallicity. We choose not to extrapolate 
beyond the grid since we found this to result in 
highly uncertain metallicities. Therefore, if any 
star falls off of the isochrone grid, either initially 
or after scattering to account for the photometric 
errors, we discard it from the final MDF. 


Since the standard approach is to assume an 
old age for most systems to which this technique is 
applied, we use isochrones for an age of 12 Gyr as 
our fiducial grid. For systems older than ~10 Gyr, 
the interpolated metallicities from RGB stars do 
not depend sensitively on the age of the isochrones 
used (Sarajedini & Jablonka 2005). Most previ¬ 
ous work with this method tends to exclude stars 
fainter than some magnitude limit in order to re¬ 
move contamination from stars on the horizontal 
branch (HB) and red clump (RC). We therefore 
limit our interpolations to magnitudes brighter 
than Mi = —1.2 mag. Finally, to obtain a singular 
estimate for the average metallicity, we construct a 
generalized histogram of the MDF (i.e. the Gaus¬ 
sian smoothed MDF using errors returned from 
the interpolation algorithm) and take the peak of 
this distribution as the representative metallicity. 
We found this to be more robust against skewing 
effects inherent in this method (see Section 3.1) 
than the mean or median. For our analysis, the 
peak of the input MDFs are calculated in a sim¬ 
ilar manner, with arbitrarily low errors assigned 
to construct the distribution. Confidence inter¬ 
vals are calculated using a standard bootstrapping 
technique, and the resulting reported errors corre¬ 
spond to 99% confidence intervals, or 3cr uncer¬ 
tainties. 


3. Results 


3.1. A Simple Stellar Population: PI 


We first synthesized the CMD of a simple pop¬ 
ulation with one age of 12 Gyr, one metallicity of 
Z = 0.001 ([M/H] = -1.28 dex), and no binaries, 
which we will refer to as PI. We apply both obser¬ 
vational profiles to this synthetic CMD, to which 
we will refer to as models lA and IB (see Section 


2.2 for a description of the observational profiles). 


The CMDs are plotted in Figure 


As was mentioned previously, studies that em¬ 
ploy interpolation on an isochrone grid typically 
limit the interpolation to stars brighter than some 
magnitude limit to prevent contamination of RC 
and HB stars on the model grid. Since a signifi¬ 
cant fraction of the AGB occupies similar magni¬ 
tudes as the bright RGB, these stars will gener¬ 
ally provide some contamination if present within 
a stellar population. Additionally, fainter stars 
have larger photometric errors, so at some point 
including fainter stars in the RGB may only act 
to add noise to interpolated MDFs. To further 
investigate how these factors affect the resultant 
MDF, we ran a series of interpolations on these 
two CMDs each time changing the faint magnitude 
limit, Mi.j, above which the interpolation was per¬ 
formed. In this way, we tested what magnitude 
range yields the most accurate and precise metal¬ 
licity estimate. We use isochrones for an age of 12 
Gyr and metallicities of Z = 0.0001, Z = 0.0003, 
Z = 0.0006, Z = 0.001, Z = 0.002, Z = 0.004, 
and Z = 0.008. The results from this set of inter¬ 
polations are presented in Table where the peak 
metallicity and its difference from the input peak 
are presented for each case. 


Overall, it is clear that the interpolations pro¬ 
duce fairly accurate values for the peak metallicity, 
given that all are within 0.1 dex of the input val¬ 
ues. We illustrate this further in Figure where 
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Table 2 

Results from varying the faint limit for the interpolation. 


Model 

Mij (mag) 

Peak [M/H] 

(dex) 

A [M/H] (dex) 

lA 

-3.5 

-1.25 

± 

0.05 

-0.03 

lA 

-3.0 

-1.24 

± 

0.06 

-0.04 

lA 

-2.5 

-1.23 

± 

0.05 

-0.05 

lA 

-2.0 

-1.22 

± 

0.06 

-0.06 

lA 

-1.5 

-1.21 

± 

0.05 

-0.07 

lA 

-1.2 

-1.20 

± 

0.05 

-0.08 

IB 

-3.5 

-1.29 

± 

0.01 

0.01 

IB 

-3.0 

-1.29 

± 

0.01 

0.01 

IB 

-2.5 

-1.29 

± 

0.01 

0.01 

IB 

-2.0 

-1.29 

± 

0.01 

0.01 

IB 

-1.5 

-1.29 

± 

0.01 

0.01 

IB 

-1.2 

-1.29 

± 

0.01 

0.01 


“^These results correspond to a simple stellar population with 
age 12 Gyr and [M/H] = -1.28 dex. 

’^Quoted uncertainties represent bootstrapped 99% confidence 
intervals. 
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Fig. 3.— Left: Difference between peak metallicity estimates from isochrone interpolation and the input metallicity 
for Models lA (circles) and IB (triangles) as a function of the adopted threshold magnitude, Mij. Right: Skew of 
interpolated MDFs as a function of the adopted threshold magnitude, Mi_/. The dashed lines represent the input 
MDF values. Top: Interpolations for Model lA. Bottom: Interpolations for Model IB. The error bars represent 99% 
confidence intervals from a standard percentile bootstrap analysis. 
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Fig. 2.— Left: The Mi, (V-I)o CMD for model lA. 
The model grid using the isochrones of [Girardi et| 
al. ( |2002 I for an age of 12 Gyr and metallicities of, 
from left to right, Z = 0.0001, 0.0003, 0.0006, 0.001, 
0.002, 0.004, and 0.008 are overplotted for reference. 
Both the AGB and RC stars are strongly blended with 
the RGB due to the photometric errors. Right: The 
CMD for model IB. In CMD IB, both the AGB and 
RC are apparent as distinct from the RGB, although 
some of these stars still contaminate the model grid of 
isochrones. 


we have plotted the difference between the input 
and output peak metallicities against Mij. 

For model IB, the interpolation proves very 
reliable owing to the smaller photometric errors 
([^ bottom left), yielding peak metallicities which 
are both accurate and precise to within 0.01 dex. 
It seems that no matter where the interpolation 
is limited, the peak metallicities in the case of 
the small error profile yield systematically lower 
metallicities. The skew of this MDF, also illus¬ 
trated in Figure]^ is consistent with skewness to¬ 
wards metal-poor values. 

While this systematic offset in metallicity is 
within the errors, we note that one possible reason 
for systematically metal-poor values could be the 
AGB stars contaminating the grid at all magni¬ 
tudes. Another possibility is that the asymmetric 
nature of the isochrone grid in the magnitude-color 
plane (i.e. the nonlinear dependence of the color of 
the RGB on [M/H]) produces this negative skew. 
Since the grid is asymmetric but the photometric 
errors are applied symmetrically, it follows that 
this effect could lead to some amount of artificial 
skewing. In order to further investigate these ef¬ 
fects, we repeated the interpolations on PI with 
the same observational profiles applied, however 
this time we excluded AGB stars from the final 
GMDs. This was achieved by applying cuts in ini¬ 
tial mass, magnitude, and color such that stars 
more evolved than stars at the TRGB were filtered 
out. Note that the cuts in magnitude and color 
were only required due to the small but nonzero 
range in the ages of the stars in this simple stellar 
population resulting from the interpolation in time 
of the SFR by lAG-STAR. We then performed 
the interpolation down to magnitudes of Miy = 
-1.2 mag. The resulting MDFs did not show any 
significant differences in peak metallicity or skew, 
forcing us to conclude that the AGB does not pro¬ 
vide a significant source of contamination for an 
old, simple stellar population. Gonsequently, the 
asymmetry of the isochrone grid remains the only 
valid explanation, and we conclude that this asym¬ 
metry does provide a significant source of artificial 
negative skew to an interpolated MDF. 

On the other hand, the model lA yields sys¬ 
tematically more metal-rich peaks than the input 
(Figure]^ top left). Additionally, this discrepancy 
grows as more faint RGB stars are included in the 
interpolation. While these differences are mostly 


6 








within the errors, its systematic nature could be 
of concern. We attribute this behavior to the in¬ 
terplay between the photometric errors and the 
isochrone grid. More specifically, the error profile 
for model lA is such that a significant number of 
stars are observed off of the isochrone grid. This 
effect becomes more prominent for fainter stars 
since they have larger photometric errors, and the 
effect is asymmetric with respect to color. That is 
to say, the grid extends further in the redward di¬ 
rection than in the blueward direction. Thus, even 
though stars are scattered symmetrically in color 
due to the photometric errors, a greater number of 
stars will ’fall off’ the model grid at the blue end. 
Simultaneously, stars that are scattered to redder 
colors are more likely to stay on the grid. This ap¬ 
pears to diminish the skew of the MDF and move 
the peak to more metal-rich values. Once again, 
the skewness for model lA (Figure top right) 
seem to corroborate this behavior. The interpo¬ 
lated MDFs of model lA appear systematically 
less skewed than those of model IB. 

Considering these effects, we now turn to deter¬ 
mining the optimal cutoff for a metallicity inter¬ 
polation on the RGB. Since we have determined 
AGB contamination to be negligible in an old 
stellar population, the fundamental competition 
governing the magnitude cutoff is between Pois¬ 
son statistics and photometric errors. On the one 
hand, excluding more stars at fainter magnitudes 
reduces the number of stars in the sample, poten¬ 
tially leading to a metallicity estimate based on a 
prohibitively small sample of stars. On the other 
hand, including stars at fainter magnitudes allows 
stars with increasingly uncertain photometry, and 
thus metallicity, into the sample. This may lead, 
for example, to the peak metallicity discrepancy 
observed for model lA. Thus, it seems that when 
the size of the photometric errors becomes com¬ 
parable to the extent of the isochrone grid, inter¬ 
polated metallicities may be subject to systematic 
errors. 

To compromise between these two competing 
effects, we examined how the number of stars in 
the interpolation sample changed with Mi j. Mov¬ 
ing the limit from Mij = -1.2 mag to Mj / = - 
1.5 mag results in losses of ^14% (lA) and ^20% 
(IB) of the stars above Mj = -1.2 mag ( note the 
difference in fractional losses is due to the com¬ 
pleteness profiles). Additionally, the interpolated 


MDFs are nearly identical in each case, with the 
only exception being a more accurate value for the 
peak metallicity in model lA (Figure]^ top left). 
Making the cutoff brighter at Mj j = -2.0 mag re¬ 
sulted in losses of 40-50% of the stars above Mj = 
-1.2 mag with no significant increase in the accu¬ 
racy of the interpolated metallicities. Under these 
considerations, we choose Mj j = -1.5 mag as our 
optimal cutoff and limit all future interpolations 
to stars brighter than this magnitude. 

An anonymous referee suggested an alternative 
method to derive the peak metallicity of a pop¬ 
ulation by constructing a fiducial RGB sequence 
from the photometry of the observed RGB. To ac¬ 
complish this, we binned all of the stars brighter 
than the optimal cutoff, Mij = -1.5 mag, into bins 
of 0.1 mag and then calculated the median (V- 
I)o color within each bin. The subsequent color 
paired with the midpoint of each magnitude bin 
served as our fiducial RGB which was input into 
the metallicity interpolation routine. We applied 
this method to both models lA and IB, and and 
the resulting differences between the input and 
output peaks in each case was 0.03 dex for model 
lA and 0.01 dex for IB. Gomparing these A[M/H] 
with those listed in Table it appears that this 
alternative method works well and recovers simi¬ 
larly accurate peaks. Flowever, since this method 
removes information about the shape of the MDF, 
we shall continue to use the original, star by star 
interpolation to construct MDFs for the remainder 
of the paper. 

Finally, we investigated how the metallicity of 
the simple stellar population affected the resulting 
metallicity interpolation. To this end, we synthe¬ 
sized a population whose properties were identi¬ 
cal to PI, but having a metallicity oi Z — 0.0076 
([M/FI] = -0.4 dex). We performed the interpola¬ 
tion on the same grid in a similar manner to PI, 
and the errors in the peak metallicity amounted to 
~0.03 dex in [M/FI]. Therefore, we conclude that 
the results of this interpolation method to be rel¬ 
atively insensitive to the bulk metallicity of the 
stellar population. 

We conclude this section with some summariz¬ 
ing remarks. For an old, simple stellar population, 
isochrone interpolation on the RGB works very 
well to retrieve the peak metal abundance, with 
metallicity errors less than 0.1 dex for the worst 
case. In cases where the photometric errors are 
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large with respect to the size of the isochrone grid, 
exemplified in model lA, the interpolated MDF 
and therefore peak metallicity are sensitive to the 
shape of the isochrone grid. For cases where the er¬ 
rors are small with respect to the model grid (IB), 
the MDF and its peak are accurately recovered 
to within 0.01 dex. We determine the best com¬ 
promise between Poisson errors and photometric 
errors of an interpolated MDF to be to limit the 
interpolation to stars brighter than M/=-1.5 mag. 
Finally, we conclude that this method is insensi¬ 
tive to the bulk metallicity of the population. 

3.2. The Effects of a Metallicity Spread: 

P2 

The limit of a simple stellar population with one 
age and one metallicity (i.e. [Fe/H]) is approx¬ 

imately valid for systems such as globular clus¬ 
ters (the effects of multiple populations and en¬ 
hancement differences in individual globular clus¬ 
ters is beyond the scope of this work). However 
more complex systems, such as dwarf galaxies, re¬ 
quire a more complex treatment. In particular, 
these systems typically show extended SFHs re¬ 
sulting in significant chemical evolution over time. 
Additionally, inherent metal abundance variations 
within a stellar population will also affect a metal¬ 
licity determination. As a simple exploration of 
the effects of a metallicity spread, we created a 
synthetic CMD whose properties are identical to 
PI (Section |3.1[ ), with the only difference being a 
spread in the stellar metallicities. To accomplish 
this, we utilize a feature of lAC-STAR which al¬ 
lows the user to provide an upper and lower bound 
for the input AMR. lAC-STAR then calculates 
the metallicity of a star formed at any instant in 
time by randomly sampling a Gaussian bounded 
by these two AMRs at that time. 

In this case, we set the lower AMR to a metal¬ 
licity oi Z = 0.0006 ([M/H] = -1.51 dex) and the 
upper AMR to a metallicity of Z = 0.002 ([M/H] 
= -0.98 dex), each constant in time, and we will 
henceforth refer to this as P2. The resulting in¬ 
put MDF is characterized by a mean metallicity 
of [M/H] = -1.20±0.15 dex, where the spread is 
calculated as one standard deviation. As with PI, 
we apply both observational profiles to this pop¬ 
ulation and refer to these two CMDs as models 
2A and 2B. We then performed the metallicity in¬ 
terpolation using the same grid from Section [XT] 


The resulting model CMDs and MDFs are plotted 
together in Figure]^ 

The peak interpolated metallicities for these 
model CMDs are [M/H] = -l.lO/lQ'Qg dex for 2A 
and [M/H] = -1.07/lo;o7 ^B, and they differ 

from their input peaks by 0 dex and -0.04 dex re¬ 
spectively. Thus, each peak is again recovered well 
within the errors. Interestingly, model 2A has its 
peak recovered more accurately than model 2B. 
Although the discrepancy is within the errors, it 
is interesting to note that the skewness of MDF 
2B is significantly more negative than for its in¬ 
put MDF. This is not true for MDF 2A. Similar 
behavior is observed for the simple stellar pop¬ 
ulation in Section 3.1 (Figure right column). 
Since the only difference between the two model 
CMDs is the observational profile applied, it seems 
that the photometric errors must be responsible 
for this. We interpret this to be related to the 
phenomenon previously observed and discussed in 
Section |3.1[ namely the interaction of the photo¬ 
metric errors with the isochrone grid. Thus, even 
though MDF 2A suffers from more noise resultant 
from the larger photometric errors, it is less neg¬ 
atively skewed than MDF 2B because more stars 
are scattered off of the blue, metal-poor end of 
the grid than the red, metal-rich end. We empha¬ 
size that this systematic effect results in recovered 
peak metallicity for MDF 2A lying near the cen¬ 
ter of the isochrone grid and is not a reflection of 
higher accuracy. This illustrates the importance 
of ensuring that the photometric errors in a CMD 
upon which metallicity interpolation will be ap¬ 
plied are small in comparison to the extent of the 
isochrone grid. In this case however, this effect is 
small. 


To add to this, the error in the peak for 2A 
is significantly smaller than for 2B, even though 
the former suffers from larger photometric errors 
than the latter. To investigate the cause of this, 
we examined the bootstrapped distributions of the 
MDF peak for both model CMDs. This revealed 
a multimodal distribution of peak metallicities for 
MDF 2B. Specifically, in addition to the true peak 
at [M/H] = -1.07 dex, there also exist less promi¬ 
nent but still significant peaks at [M/H]~-1.22 dex 
and [M/H]~-1.03 dex. Careful examination of the 
generalized histogram of MDF 2B (Figure]^ bot¬ 
tom right) reveals these secondary peaks. These 
secondary peaks act to add uncertainty to the pri- 




Fig. 4.— Left: The Mi, (V-I)o CMD for models 2A (top) and 2B (bottom). The model grid using the isochrones 
of Girardi et al. (20021 for an age of 12 Gyr and metallicities of Z = 0.0001, 0.0003, 0.0006, 0.001, 0.002, 0.004, and 
0.008 are overplotted for reference. Right: The interpolated binned (points) and Gaussian smoothed (dashed line) 
MDFs for each corresponding model CMD. The error bars in the binned MDFs are la Poisson uncertainties. The 
solid bar plots represent the input MDFs. 
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mary peak via our bootstrapping technique. It is 
therefore likely that we are over-estimating the er¬ 
rors in this case by using the 3 ct errors. On the 
other hand, the larger photometric errors in CMD 
2A smear these peaks out beyond recognition such 
that they do not present themselves in the boot¬ 
strapped distribution of the MDF peak. Taking 
into consideration all of the effects discussed in 
this subsection, we conclude that the peak inter¬ 
polated metallicity is relatively insensitive to the 
effects of an inherent metallicity spread of order 
cr ~ 0.15 dex. 

We conclude this subsection by noting that 
overall, the shape of the interpolated MDF 2B (|^ 
bottom left, points) closely matches that of the in¬ 
put MDF (|^ bottom left, solid line). The same is 
not true for MDF 2A Q. We note that the input 
MDFs for the two model CMDs differ slightly due 
to Poisson statistical effects introduced during the 
application of the completeness profile. 


3.3. The Effects of an Age Spread: P3 


We now turn to exploring the effects of a spread 
in the stellar ages. For this we synthesized a syn¬ 
thetic CMD whose properties are similar to the 
simple stellar population, PI, but whose SFH con¬ 
sists of a Gaussian with a mean age of 12 Gyr and 
a Icr width of 2 Gyr. We designate this popula¬ 
tion as P3, and again apply the two observational 
profiles creating two models GMDs 3A and 3B. 
We performed the interpolation as before, and the 
results are illustrated in Figure 


The peak interpolated metallicities in this case 
are [M/H] = -1.20jlQ;o5 ^A and [M/H] = 

-1.30i11q;q 2 fo'' 3B. We note that these results 
are similar to that for the simple stellar population 
PI (see Section 3.11. That is, the peak metallicity 
is well-recovered in the case of smaller photometric 
errors, albeit slightly metal-poor. The only effect 
of the age spread in this case is to increase the 
scatter from ^0.01 dex to a few hundredths of a 
dex. On the other hand, the peak for 3A with its 
larger photometric errors is again more metal-rich 
than the input metallicity. This seems to be the 
result of the interplay of the photometric errors 
with the isochrone grid. In this case, the uncer¬ 
tainty in the peak remains the same as for a simple 
stellar population, indicating that the photomet¬ 
ric errors wash out any signal of the age spread in 


this case. We therefore conclude that for an old 
age, the effects of an age spread of order cr ~ 2 
Gyr do little to affect a metallicity interpolation. 
The peak is recovered with similar accuracy to a 
single-age population. The only observed effect is 
an increase in the error of the peak for the small 
error case, resulting in uncertainties still well be¬ 
low 0.1 dex in both cases. 

3.4. The Effects of an Age and Metallicity 
spread: P4 

We now add complexity to the synthetic GMD 
by inputting spreads in both metallicity and age. 
For this case, we essentially combine the AMR 
from P2 and the SFH from P3 and create an old 
population with mean age 12 Gyr, age spread of 
2 Gyr, mean metallicity of [M/H]= —1.20 dex, 
and metallicity spread of 0.15 dex. This popula¬ 
tion is designated P4, and applying the observa¬ 
tional profiles as before we have two model CMDs 
4A and 4B. Performing the interpolations in the 
same way as with the previous models, we con¬ 
struct the MDFs along with their corresponding 
CMDs in Figure 

The peak interpolated metallicities for mod¬ 
els 4A and 4B are [M/H] = -1.04/lg Qg dex and 
[M/H] = -1.07jlg;|!)5 dex respectively. The peaks 
for both input MDFs are near [M/H] = -1.00 dex, 
so both peaks are recovered within their respec¬ 
tive confidence intervals. One may note that the 
quoted uncertainties for these values are similar 
to those of models 2A and 2B, indicating that the 
metallicity spread dominates over the age spread 
in this case. Again, here it appears that we may 
be overestimating the errors in the case with the 
smaller photometric errors. The same mechanism 
discussed for model 2B in Section [32] seems to be 
at work for model 4B, namely the asymmetry of 
the isochrones in the CMDs and how this affects 
the metallicity interpolation. Aside from this, the 
peak metallicities are well-recovered in this case 
with both an age and metallicity spread. 

3.5. The Effects of an Age Error 

In the absence of an independent measure of 
the age of a stellar system, it is also interesting 
to investigate how a systematic error in the as¬ 
sumed age will affect a metallicity determination 
using the RGB. As was discussed above, for an 
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Fig. 5.— Left: The Mi, (V-I)o CMD for models 3A (top) and 3B (bottom). The model grid using the isochrones 
of Girardi et al. (20021 for an age of 12 Gyr and metallicities of Z = 0.0001, 0.0003, 0.0006, 0.001, 0.002, 0.004, and 
0.008 are overplotted for reference. Right: The interpolated binned (points) and Gaussian smoothed (dashed line) 
MDFs for each corresponding model CMD. The error bars in the binned MDFs are la Poisson uncertainties. The 
solid bar plots represent the input MDFs. 
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Fig. 6.— Left: The Mi, (V-I)o CMD for models 4A (top) and 4B (bottom). The model grid using the isochrones 
of Girardi et al. (20021 for an age of 12 Gyr and metallicities of Z = 0.0001, 0.0003, 0.0006, 0.001, 0.002, 0.004, and 
0.008 are overplotted for reference. Right: The interpolated binned (points) and Gaussian smoothed (dashed line) 
MDFs for each corresponding model CMD. The error bars in the binned MDFs are la Poisson uncertainties. The 
solid bar plots represent the input MDFs. 
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Fig. 7.— The model CMDs used in investigating the effects of an age error on a metallicity interpolation. 
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old population (>10 Gyr) the RGB is relatively 
insensitive to the exact age. However, for younger 
populations this does not hold. To explore this 
further, we generated a series of synthetic CMDs 
each with the same properties as P4, except chang¬ 
ing the mean age of the stars by 1 Gyr from 11 
Gyr down to 3 Gyr. We then applied the obser¬ 
vational profiles and interpolated in metallicity on 
the 12 Gyr isochrone grid to observe how the peak 
metallicity varies with a fixed adopted age of 12 
Gyr. Each of these model CMDs is displayed in 
Figure 

Figure shows how the metallicity discrepancy 
changes as a function of mean population age for 
model CMDs using the observational profiles A 
(circles) and B (triangles). Here, A[M/H] is in 
the sense [M/H]out -[M/H]in. Figure shows a 
clear trend in which populations younger than the 
isochrones in the model grid produce systemati¬ 
cally more metal-poor metallicities. In accordance 
with previous studies, this effect is small (^0.1 
dex) in populations older than 10 Gyr but may 
produce an error of nearly half a dex for consider¬ 
ably younger populations. Linear fits to the data 
resulted in the parameters presented in Table 
The relations are incompatible with each other 
considering the quoted uncertainties. It seems 
that CMDs with smaller photometric errors along 
the RGB are more sensitive to age differences in 
their interpolated MDFs. Part of the decreased 
sensitivity for profile A may be the coupling be¬ 
tween the isochrone grid and larger photometric 
errors discussed in Section [tTI 

Of particular note for the small error profile 
(triangles) is the discontinuous jump in metallic¬ 
ity error from 6 to 5 Gyr. It seems that this age 
range defines where the stars present on the inter¬ 
polation grid become increasingly dependent on 
the age of the system. 

On the other hand, for the case with larger 
errors (circles) we observe systematically more 
accurate peak metallicities. While this is cer¬ 
tainly counter-intuitive, this is almost certainly 
a result of the coupling between the photometric 
errors and the adopted isochrone grid. As dis¬ 
cussed in Section |3.H the large photometric er¬ 
rors in profile A couple with the asymmetry of 
the isochrone grid to scatter more stars off of the 
metal-poor end compared with the metal-rich end 
of the grid. Thus, the seemingly more accurate 
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0.0 


- 0.1 

S’ 

i; -0-2 

<1 

-0.3 

-0.4 
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Fig. 8.— Peak metallicity discrepancy as a function 
of the mean age of a synthetic CMD. Circles represent 
CMDs modeled with observational profile A, while tri¬ 
angles are CMDs modeled with profile B. The dashed- 
dot and dotted lines show the linear fits for the two 
profiles A and B, respectively. 
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Table 3 

Linear fit parameters for metallicity error versus age offset experiments. 


Profile 

Slope‘s (dex/Gyr) 

Zero-point®' (dex) 

RMS scatter (dex) 

A 

0.020T0.007 

-0.235±0.055 

0.049 

B 

0.039±0.010 

-0.501±0.079 

0.069 


’^Quoted uncertainties represent the square roots of the diagonal ele¬ 
ments from the covariance matrix returned from the Htting routine. 


values for CMDs obtained applying profile A is 
more likely an artifact of the interpolation proce¬ 
dure rather than a reflection of the true accuracy. 

3.6. Multiple Gaussian bursts: P5 

We tested the effects of multiple epochs of star 
formation in a system by simulating a stellar popu¬ 
lation with 3 distinct star formation episodes. We 
achieved this by inputting a SFH into lAC-STAR 
which consisted of three Gaussian bursts of star 
formation at ages of 12, 8, and 4 Gyr. All have 
Icr widths of 1 Gyr, and the AMR is monotoni- 
cally increasing between the peaks of star forma¬ 
tion. Both the SFH and AMR for this popula¬ 
tion, which we call P5, are illustrated in Figure]^ 
Applying the observational profiles and interpolat¬ 
ing on the same isochrones as in previous sections 
yields CMDs 5A and 5B. These CMDs and cor¬ 
responding interpolated MDFs are illustrated in 
Figure [T0| 

The peak of the input metallicity lies at [M/H] 
= -0.98 dex. Since the most recent star forma¬ 
tion episode at 4 Gyr dominates the number of 
stars currently present on the RGB, the MDF is 
peaked at this metal-rich end. However, the inter¬ 
polated MDFs for both CMDs 5A and 5B are sig¬ 
nificantly more metal-poor than this with peaks at 
[M/H]= -1.181°;°^ dex and [M/H]= -1.22^1^ 
dex, respectively. It therefore seems that the most 
recent episode of star formation is being mistaken 
for a more-metal poor population owing to the old 
isochrone grid being used. This is in accordance 
with the results presented in Section |3.5[ where 
we showed that for systems younger than ~6 Gyr, 
the shape of the RGB becomes heavily dependent 
on the age, and relatively large errors in the peak 
metallicity are expected. 


It is interesting to note that in the case of the 

lower 


smaller error profile (CMD 5B, Figure 10 


right), two distinct peaks in the MDF are recov¬ 
ered through the interpolation. These peaks corre¬ 
spond roughly to the metallicity of the system dur¬ 
ing the peaks of the first two episodes of star for¬ 
mation ([M/H]^-1.5 dex, 12 Gyr burst; [M/H]^- 
1.3 dex, 8 Gyr burst; see Figure [^. These peaks 
are well-recovered in the interpolated MDF, high¬ 
lighting the high fidelity of the RGB for metal¬ 
licity estimates of old stars. On the other hand, 
the metallicity of the system during the youngest 
starburst ([M/H]~ -1 dex, 4 Gyr; see Figure 
completely lacks a peak in the interpolated MDF, 
further illustrating the ineffectiveness of the RGB 
metallicity estimation technique for stellar popu¬ 
lations younger than 6 Gyr absent of an accurate 
age estimate. 


3.7. The Effects of a Recent Star Forma¬ 
tion Episode 

We now turn to investigating the effect of a 
young starburst in more detail. To accomplish 
this, we synthesized a series of stellar populations 
each with two Gaussian star bursts, one at 12 Gyr, 
another at 4 Gyr, and both with Icr widths of 1 
Gyr. In each case of this series, the intensity of 
the young star formation episode was scaled down 
by a multiplicative factor such that its intensity 
ranged in 0.25-1.0 times the intensity of the old 
star burst. The AMR for all populations in this 
series are identical to P5. The same procedures for 
applying the observational profiles and performing 
the metallicity interpolation were then performed 
for each of these GMDs. Figure pi] shows the error 
in the peak metallicity of the interpolated MDF as 
a function of young starburst intensity. 
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Fig. 10.— Left: The Mi, (V-I)o CMD for models 5A {top) and 5B (bottom). The model grid using the isochrones 
of Girardi et al. (20021 for an age of 12 Gyr and metallicities of Z = 0.0001, 0.0003, 0.0006, 0.001, 0.002, 0.004, and 
0.008 are overplotted for reference. Right: The interpolated binned (points) and Gaussian smoothed (dashed line) 
MDFs for each corresponding model CMD. The error bars in the binned MDFs are la Poisson uncertainties. The 
solid bar plots represent the input MDFs. 
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Relative intensity of young starburst 


Fig. 11.— Peak metallicity error versus intensity of a 
young starburst (4 Gyr) relative to the old star burst 
(12 Gyr). CMDs produced with observational profile 
A are circles, while the triangles show CMDs produced 
with observational profile B. Error bars represent 99% 


12 10 8 6 4 2 Oconfidence intervals returned from bootstrapping the 


Look-back time (Gyr) 


peak of the interpolated MDF. 


Fig. 9.— Top: The input SFH for P5 consisting of 
3 Gaussian bursts of star formation. Note that lAC- 
STAR cannot take more than 20 nodes in time for 
the input SFH, hence the coarse resolution here. lAC- 
STAR does however interpolate amongst these nodes 
in time to achieve greater time resolution. Bottom: 
The input AMR for P5. 
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Inspection of Figure 11 reveals that for young 
starbursts <0.3 as strong as the old star burst, the 
peak of the MDF is well-recovered by interpolat¬ 
ing on a grid of old isochrones. Once this young 
burst exceeds this threshold intensity, the young 
stars begin to dominate the MDF and the peak in¬ 
terpolated metallicity is in signihcant error when 
compared to the peak of the true MDF. Taking 
these results with those of the previous subsec¬ 
tion, we conclude that a metallicity determination 
from the RGB of a stellar population that has ex¬ 
perienced star formation within the last ^6 Gyr 
that was at least 30% as intense as star formation 
in older epochs will produce a metallicities which 
are systematically in error by up to 0.5 dex. 


3.8. The Effects of Binary Stars 

Another component of stellar systems that may 
act to alter an interpolated MDF is binarity. Stud¬ 
ies of the SFH of dwarf galaxies in the Local Group 
typically take the binary fraction as a free parame¬ 
ter in their analysis, and they report values around 


40-60% ( 

Hidalgo et al. 

2009 Monelli et al. 

2010 

Hidalgo et al.|2011 

|. We choose to test the affects 


of binary systems on the metallicity interpolation 
by introducing a binary population to PI. This 
population contains binary stars in 40% of the sys¬ 
tems with a minimum mass ratio of q = 0.5, with 
the mass of the secondary chosen from a flat mass 
distribution. Subsequent interpolation as with the 
other synthetic GMDs did not reveal noticeably 
different MDFs from PI, forcing us to conclude 
that stellar populations with binary characteristics 
similar to those used here should not significantly 
affect a metallicity interpolation on the RGB. 


4. Conclusions 

We have explored a variety of different stellar 
population effects on the MDF constructed from 
the RGB stars in synthetic GMDs. In particu¬ 
lar, we examine how realistic photometric error 
and completeness effects coupled with a variety 
of SFHs and AMRs alter the metallicities derived 
from an interpolation in Mj, (V-I)o, and Z on the 
RGB. We summarize our findings for a variety of 
simulated model GMDs as follows: 


1. For an old, simple stellar population, the 
RGB proves to be a very reliable metallicity 


indicator. In the worst case, peak metallic¬ 
ity values are recovered within less than 0.1 
dex for photometry with errors that are rela¬ 
tively large compared with the extent of the 
isochrone grid. In the best case, peak metal¬ 
licities are recovered to within 0.01 dex for 
more precise photometry. We determine the 
optimal magnitude range of RGB stars over 
which a metallicity interpolation should be 
performed to be restricted to stars between 
the RGB tip and Mj = —1.5 mag for the 
range of the photometric errors used in this 
exploration. Inherent spreads in the stellar 
metallicities and ages do not have signihcant 
effects on the modular interpolated metallic¬ 
ity. 

2. It should be noted that the adopted isochrone 
grid can introduce signihcant systematic ef¬ 
fects in the resultant MDF. Signihcant se¬ 
lection biases may present themselves in the 
subsequent interpolated MDF resulting from 
asymmetric scattering of the stars off of the 
model grid. 

3. A stellar population signihcantly younger 
than 10 Gyr will yield systematically more 
metal-poor metallicities assuming an old age 
for the interpolating isochrones. The re¬ 
sulting discrepancy in peak metallicity as a 
function of age can be described with a lin¬ 
ear relation. The slope of this relationship 
depends on the photometric error prohle of 
the GMD, indicating that more precise pho¬ 
tometry is more sensitive to this age error 
than photometry with larger errors. Thus, 
care should be taken to obtain an accurate 
age estimate when estimating metallicities 
from the RGB using precise photometry, es¬ 
pecially if a signihcant possibility exists that 
the population is younger than 10 Gyr. 

4. Using a grid of old isochrones, a star for¬ 
mation episode occurring less than ~6 Gyr 
ago produces a signihcant, systematic error 
in the interpolated MDF. This can result in 
metallicities that are erroneously poor by up 
to 0.5 dex. Along the same vein, if a recent 
starburst is >30% as intense as older star 
formation events with comparable durations 
in time, the interpolated MDF will yield a 
signihcantly erroneous peak. 
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We close by noting that the RGB of old stellar 
populations has proven itself as a reliable metallic- 
ity estimator in our study. For reasonably simple 
SFHs, the RGB provides a very accurate estimate 
of the average metallicity of a stellar population. 
However, if the system in question has experienced 
any significant star formation within the past few 
Gyr, metallicities derived from the RGB may be 
in error by up to 0.5 dex. Therefore, when ap¬ 
plied to stellar systems with potentially complex 
stellar populations, more care needs to be taken if 
using the RGB to estimate the metallicity of the 
system. In particular, a reasonably accurate esti¬ 
mate of the age of the system should be obtained 
through an independent, reliable estimator. Addi¬ 
tionally, one should be reasonably certain that no 
significant star formation has occurred within the 
past few Gyr in the system in question, otherwise 
there is a high risk of significant systematic errors 
in the derived metallicities. 
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